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Abstract 

With the aim of exploring a massive model corresponding to the perturbation 
of the conformal model jlj the nonlinear integral equation for a quantum system 
consisting of left and right KdV equations coupled on the cylinder is derived from 
an integrable lattice field theory. The eigenvalues of the energy and of the transfer 
matrix (and of all the other local integrals of motion) are expressed in terms of 
the corresponding solutions of the nonlinear integral equation. The analytic and 
asymptotic behaviours of the transfer matrix are studied and given. 
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1 Introduction 



The computation of the eigenvalues of the conserved charges - in particular the local 
ones - is a very important issue in the study of 1 + 1 dimensional integrable models. For 
instance, the local integrals of motion have recently revealed to enter the form of off- 
critical null- vectors [2]. Among the various methods of computing, a possible one uses 
a discretisation of the space and the diagonalisation of the discrete transfer matrix by 
means of a Bethe Ansatz technique. The problem of recovering the continuous limit of 
such eigenvalues can be tackled using nonlinear integral equations, method previously 
applied to genuine lattice models like the XXZ chain and the inhomogeneous six- vertex 
model [31113]. 

In a previous paper ^ we coupled on the lattice the Lax operators of two - left 
and right - quantum (m)KdV equations, proposing two slightly different monodromy 
matrices as possible descriptions of perturbed minimal conformal field theories. The 
Bethe Ansatz diagonalization of the discrete transfer matrices was completely achieved, 
showing that the problem is somehow equivalent to a twisted spin —1/2 XXZ chain 
with alternating dishomogeneities. The twist was called dynamical, because it is au- 
tomatically generated by the theory (and depends on the number of Bethe roots), not 
introduced as an external ad hoc parameter. We concluded that this lattice model de- 
scribes, in some sense, a dynamically twisted version of the lattice sine-Gordon theory 
[7j and, in this respect, we conjectured and partially argued that in the continuous limit 
it gives a description of perturbed minimal conformal models. 

In this letter we address the problem of performing the continuous limit of the 
eigenvalues of the aforementioned transfer matrices |W , by finding the nonlinear integral 
equation obeyed by the counting function. In particular, exact expressions of all the 
local integrals of motion are given as nonlinear functionals of the counting function. 
This problem is different from that solved in jH |3] , since these papers used a lattice 
approach to the twisted sine-Gordon model based on an inhomogeneous spin +1/2 chain 
(with a suitably added twist). However, this lattice discretisation does not furnish a 
transfer matrix for the sine-Gordon model, whereas our construction is endowed with 
an off-critical transfer matrix, which generates all the integrals of motion, in particular 
the local ones. On the other hand, our lattice model satisfies a braided version of the 
Yang-Baxter algebra [H]. 

2 The nonlinear integral equation for the vacuum 
sector 

In this paper we continue the work done in [ij, extending it to the massive case in an 
unifying view of conformal and off-critical theory. In fact, we start by recalling the 
introduction of the mass 0, realised through coupling two KdV theories via dishomo- 
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geneities. More precisely, we considered on a lattice of length R, spacing A and number 
of sites = -R/A, G 4N, the following monodromy matrices 



N/i 



M(a) = e-"^ JjL4.(a + e)L4i_i(a + e)L4i_2(a-e)L4i_3(a-e), (2.1) 



1=1 

N/4 



M'{a) = e"^ JjL4i(a-e)L4i-i(a-e)L4i-2(a + e)Z4.-3(a + e) , (2.2) 



i=l 



with coupling constant < < 1 and dishomogeneity real. In these formula left 
Lax operators Lm{<y) and right Lax operators Lm{c() are given by 



-iVm 



Lmia) 



(2.3) 



where are the discretised quantum counterparts of the mKdV variables v (left case) 
and V (right case) and satisfy the nonultralocal commutation relations 



\v+ , V- 



(•^m— l,n 2(5ni,n ~l~ ^m,n—l) 



(2.4) 
(2.5) 



and the periodicity conditions V^+jy = V^. In fact, we have proved in P that, as a con- 
sequence of ()2.4t 12. 5p . monodromy matrices ()2.1|) and ()2.2p satisfy braided Yang-Baxter 
equations [S]. Moreover, the diagonalisation of the corresponding transfer matrices [H] 
is based on the Bethe equations 



n 



sinh(as — + ^7r/5^ 
^ sinh(as — — ztt/^^ 



sinh ( a, + e - sinh f a, - 6 - ^ ' 



sinh ( a, + e + ^ ) sinh (a, - 6 + ^ 



(2.6) 

in which / is the number of Bethe roots and e is equal to +1 for ()2.H) and to —1 for 
()2.2|) . Correspondingly, the eigenvalues on Bethe states of the transfer matrices are 



A^(a) = A+(a) + A^(a) 



(2.7) 



where 



and 



sinh(a — a,. ± ivr/?^) . 
^ smh(a — ar) 



eiv 



4 sinh 6 - q; =F — r— ) sinh | 6 + a ± 



N/4: 



(2.8) 



(2.9) 
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Comparison with equation (2.18) of jH] shows that, after the identifications 2cosh2G = 
1/ S, = As — m /2 and the substitution /?^/8vr, Bethe equations ()2.6|) coincide 

with the Bethe equations for the lattice sine-Gordon model jTj twisted by the dynamical 
factor e-^2i7r/32«^ 

Now we want to summarise all the Bethe equations ()2.(ij) into a single nonlinear 
integral equation. As a preliminary, we define the function, analytic in the strip |Imx| < 
min{C,vr-C}, < C < vr, 

^, sinh« + x) 

0x,C^Un— 2.10) 

smh(2(^ — x) 

which allows us to define the counting function 



7 ( 

Zn[x) = — 



(2.11) 

In terms of the counting function the equations ()2.fjj] assume the form 

ZN{as) = 7r(2/s - / - 1) , /. G N . (2.12) 

Between the various solutions to the Bethe equations we single out the vacuum, i.e. the 
eigenstate with the minimum energy, for which are real and Z7v(«s) are equal to all 
the numbers of the form ()2.12|) between Z]\i{—oo) and Zjv(+oo). From ()2.11|) it follows 
that these two extremes are 

Zn{±oo) = ±^(27r - 27rp^) ± /(vr - 27ip^) + 27rP^el . (2.13) 

The vacuum minimises the energy - which we will define in next section - when 

l/3<p^<2/3, (2.14) 

and the "twist" 27r/5^e/ is small enough mod. 27r. Considering in the region ()2.14j) . 
we now want to write an equation for the counting function associated to the vacuum 
in the continuous limit: +oo, R fixed. We firstly remark that ()2.13|) imply that 

in the continuous limit also / +oo and therefore, even if we define 

P'^ = p/p' , p < p' coprimes , (2.15) 

the last term in ()2.11|) is divergent. Therefore, we stick to the different subsequences 

l = np' + K , 0<K<p'-l, (2.16) 

where n +oo with k fixed. Discarding in the definition of Zn{x) multiples of 27r we 
get: 

I 

Zn{x) = - [0 (^x + e, + (^a; - e, ^Z^^)] + J] (a; - a,, vr/J^) +27r£a; , (2.17) 



r=l 
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where uj is finite in the continuous hmit (the double braces denote the fractional part): 

Now, using standard techniques (see |4^), one can rewrite a sum over the vacuum Bethe 
roots of a function with no poles on a strip around the real axis as follows: ^ 



Y,fN{ar) = - ^/;(x)Z^(x) + 2 / ^/^(x)Imln[l + e^^-(^+^°)] . (2.19) 



Applying ()2.19j) to the sum contained in ()2.17|1 we get 



(^x + e, + (^x - 0, ^Z?')] + 27Teuj + (2.20) 

t/ — OO — CXD 



Let us define for sake of conciseness 

K(x) = ^<P' {x, 7r/32) , $(x) = (x + e, f + (x - 6, f/^^) , 

L^(x) = Imln [1 + e^^^(^+^o)] . (2.21) 

In these notations the Fourier transform of relation ()2.20p is 

ZNik) = V 2 V^Liv^ +47r2^c^ M — , 2.22 

4 l-ir(A;) l-K{k) 1-K{k) 

where we have adopted the convention f{k) = dxe~^^^ f (x) about the Fourier 
transform. Explicitly, the Fourier transforms of the first two functions in ()2.21|) are 
respectively 

. , sinh (f - 7r/32) sinh (^ - ^Z^^) 

K(k) = \' , ^ , $ A; = \' ^ cos ke . 2.23 

^ ^ sinhp ' ^ ' ik sinhffc ^ ' 

Reparametrising the coupling constant introducing the kernel G{9) and the renor- 
malised twist u), 

7-vr(l-/5^), G[k)^- ^('^ = sinhfc(|-7) . 
' ^ l-ir(A:) 2coshfc|sinhA;(f -|) ' /32 ' ^ ^ 

equation ()2.22j) takes the compact form 

Z^ik) = —- — + 2G(A;)Ljv(A;) + 2neC^5{k) . (2.25) 



^We are fixing I to be even in the following (i.e. e'^"'"^^ = —1), as the case I odd can be treated 
along the same lines and in the end does not add new information. 
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Now we want to Fourier anti-transform and then to perform the continuous hmit: 
therefore, it is necessary to work out in such a hmit the behaviour of the anti-transform 
of the first term in the r.h.s. 

Fn{x) = 

J-oo l-K[k) J 

ik{x+e) _L ^ik{x-e)-i sinhfc (| - |/?^) 

^ sinhA;f/32coshA;(f - f/^^) ' ' 

This driving term may be evaluated through the residue theorem and it diverges in 
the continuous limit, unless 9 depends on in a peculiar way. We sketch here the 
evaluation, since it differs from the previous ones IH E] (even from the conformal case 
PP). All the poles of the integrand lie on the imaginary axis and, when 1/3 < jS"^ < 2/3, 
the poles with smallest modulus are k = 0, ±z^. If in the continuous limit G — > +oo, we 
evaluate for fixed x the integral in ()2.26|) containing e'^^C^+Q) by choosing a semi-circle 
contour in the upper fc-complex half-plane with a small semicircle avoiding the pole at 
A; = 0. Since the contributions from the other poles are damped by an exponential of G, 
the leading term is given by the residue at the pole with the smallest nonzero modulus, 
i.e. k = iii'/'j, minus the contribution of the semicircle around k = 0, whose value is 
linear in A^: ^^§^"2^ Analogously, in the evaluation of the integral containing e*'^^^"®) 
we close the contour - containing a semicircle avoiding = - in the lower fc-complex 
half-plane. Now, the leading contribution is given by the pole a.t k = —iTr/'j, minus 
the contribution of the semicircle around k = 0, whose value ( — I cancels the 



previous one. Adding the two residue contributions we obtain (up to terms o{N^) — 
when A^ +00): 

FN{x)=Ne~^^- ^^ + o(A^°) =mi?sinh-x + o(A^°), (2.27) 

sm 2(1-/32) 7 

where we have supposed a logarithmic divergence of G, via a constant mass parameter 
m (in such a way that niR is dimensionless) : 



G = --In 



sm 



A^ 2(1 -/?2) 



(2.28) 



With this result in mind it is easy to write an equation satisfied by the counting function 
in the continuous limit, Z{x), from relation (j2.25p : 

Z(x) = mR sinh -x + 2 dy G(x - y)lm In [l + e*^(^+*°)l + eu . (2.29) 

7 J-oo 

This nonlinear integral equation describes the tu-vacuum of the theory at the different 
values of Co. 
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Remark 1 Equation ()2.29p has been found for 1/3 < (3"^ < 2/3 and small u. 
However, it can be considered without problems in the whole region < < 1 and 
u} > 0, defining by analytical continuation a state of the continuous theory which 
minimises energy. 

Remark 2 Equation ()2.29|1 with u = would also describe the vacuum of the 
continuous theory from the lattice sine-Gordon model [31^], with lagrangian 

^ = 2 + ^ ' ^2-2°) 

where mo = 2\/2 ^ (cosh 20)"^/^ is a bare mass^. Hence, there is a sort of inter- 
change, given by the first of ()2.24|) . between the coupling constant (3 in the lagrangian 
and 7 in the equation (|2.29|) . 



3 Local charges and transfer matrix 

The continuous limit of a sum like (|2.19p can be expressed, once (j2.29|) has been put in, 
in terms of only the counting function Z{x) (and the continuous limit of /Ar(a;), f{x)): 

hm y2f^(ar) = ft+ / — J/(x)Imln[l + e^^(^+^°)] . (3.1) 

NA=Rr=l 



We have rearranged the r.h.s. in order to separate the bulk contribution 



fb = ^ lim 



— [ —f'j^{x)(mRsinh—x + eu 
J-oo 2tt \ 7 



(3.2) 



from the finite size corrections depending on the linear functional 

dk 



J fix) = lim 



—e'^-tkfN{k)J{k), J{k) 
Zn 



sinh fc| 



2 cosh sinh ^ (f — ^ ) 



(3.3) 



Formula (j3.ip allows to calculate exactly state depending quantities of a Quantum Field 
Theory on cylinder. In particular this formula applies effectively to the calculation of 
the energy of the tD-vacuum. We build up the hamiltonian operator in tight analogy 
with the same construction in lattice sine-Gordon model 



H 



1 



4A sin7r/52 cosh 26 



e-^'^—G-Ja 



d 



da 



,2a 



d 



d 



a=e+ 



(3.4) 



^From (|2.28|l we have the renormahsed mass parameter m = 
continuous hmit. 



mo ( R\l^ 
4 \N I 



sin 2(7-/3^) the 
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where we have defined = In [p^(a)~^A7v(a)] • Using ()2. 7112.91 13.4|) . the eigen- 

i 

values of H for a general Bethe state read as a sum on Bethe roots, Yl h{ar, N, R), 

r=l 

where 



h{x, N, R) 



4A sin 7r/?2 cosh 26 I gj^h (^x + & + ^) ^ sinh (^a; + 9 + ^ 



■se-x- 



sinh (x-e + ^) sinh (^x - 9 + ^) sinh (^x + 9 



rv(3-5) 



+ 



sinh ( a; + 9 — 



+ 



sinh ( a; — 9 — 



sinh f X — 9 — 



i7r/32 



When 1/3 < 13^ < 2/3, extending the results of |IU] to u; ^ 0, we are able to show that 
the cD-vacuum is indeed the eigenstate with the lowest eigenvalue. Moreover, when the 
twist a) = this is the energy of the lattice sine-Gordon model [7j. 
After Fourier transforming ^T] 



ikh{k,N, R) 



20 



vrfcsinh A;|/3^ 



sin [m - TT/^^) e-^l^^-l^' 



A cosh 29 sin 7r/5^ sinhA;| 
- sin {kQ + 7r/52) e'^f-'^^^l , 

we are left with the evaluation of Jyjyx) ()3.3p in the continuous limit 



(3.6) 



lim 



20 



N^oo [2_Rcosh29sin7r/92 
- sin (A;9 + 7r/32) e'^t-'^"^' 



dk 



ikx 



Tvk 



2t[ coshfc^ 



sin (A;9 - vr/?') e-^'i+'="^' 



(3.7) 



Again, because of relation ()2.28|) we can calculate exactly this limit by using the residue 
method (closing the contour of integration in the upper or in the lower /c-complex half- 
plane). Collecting the finite contributions we obtain a simple behaviour 



sm 



Jh(x) 



-4m- 



TT 1-2(3* 



2 1-/32 



2(1-/32) ^ 

1 ^2^2^ — smh-x. 

1 — p^)"^ siuTrp"' 7 



(3.^ 



Hence, the finite size correction E ~ Ei, to the energy of the cj-vacuum is 

7r/32 



sm 



E-Eh 



-Am- 



7T 1-213* 
2 1-/32 



sm 



2(1-/32) 



— Sinh -xlmln h + e^^(-+^o)l 



TT 



(3.9) 



As far as the ci;-vacuum state is concerned, this formula has been derived within a 
twisted lattice sine-Gordon model and hence proves the quantum equivalence with the 
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approaches to the sine-Gordon model based on spin +1/2 XXZ chain 01 However, 
a different renormahsation of the mass occurs in ()3.9|) . 

We now want to compute the continuous hmit of the eigenvalues of the transfer 
matrix ()2.7|) proposed in 6j as description of perturbed conformal field theories. To 



make formula (jH.lll useful, we shall define F 



N 



a 



InA^(a). Let us firstly concentrate 



on F+(a) 



FJ(a) 



sinh(a — + inP^ 



N 



ln2sinh i 6 



a 



sinh(a — a 

in (3^' 



+ 



+ In 2 sinh 6 + a + 



171(3' 



20 



. (3.10) 



Now we restrict ourselves to the u)-vacua, although generalisations to excited states 
can be worked out along the developments of [S]. The last addendum in ()3.10p gives a 
contribution whose continuous limit is zero if < < 1/2 and infinity if 1/2 < < 1. 
Therefore, we shall regularise F"'"(a) for 1/2 < < 1 defining it as the analytic 
continuation of F"'"(a) with < < 1/2. With this regularisation in mind we obtain 
from dnH) 



F+(a) = F+(a) + 



dx 



dk i 



27r 



e'^'^ikf+ik, a)J{k) Imln [l + e^^(^'+*°)] - eiixoo , 

(3.11) 



where F^(a;) is the bulk contribution, J{k) is given by ()3.3|) and 

/+(a;,a) = In ■ 



sinh {a — X + ivr/?^) 



(3.12) 



sinh(a — x) 

The twist uj comes out as in [1 . Since F+(a) ()3.10p has the periodicity property 
F+(a + ivr) = F+(a;), it is sufficient to study F+(a) in the strip 7 — vr < Ima < 7. 
Comparing p.lip with relation (5.6) of [1^, we remark that the second term in the r.h.s. 
is formally the same. Therefore, we can use the results of that paper: from formulae 
(5.16, 5.24) of Ij we obtain that, if < Ima < 7, 

"+~ dx 1 



7 sinh ^ (x 



a 



-Imln [l + e'^("+'°)] -dvrcu, (3.13) 



and that, if 7 — vr < Ima < 0, 



a) — emu + 



ik\ x—a—- 




i sinh k^ 



2 cosh sinh ^ (f — f) 
On the other hand, from definitions ()3.2|) and ()3.12|) it follows that 
dx sinh z7r/3^ 



(3.14) 

Imln[l + e'^("+'°)] . 



F+l 



a) 



27r cosh (a — a; — y) cosh (a — x — y + ivr/J^) 



TT 



mR sinh —x + eu . 
7 

(3.15) 
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This integral is finite only for 0</5^<l/2. IfO< Imo; < 7 its value is 

,2 



-mR ctg — cosh —a + lenu 
27 7 



(3.16) 



whereas, if 7 — tt < Ima < it equals 



F+(a) = -mi?-^ 
31112:^ 



cosh 



"tt / 


.7r\ 


— a - 




_7 V 


2 



+ isT^uj — ieu 



(3.17) 



Let us now consider F (a) which reads from ()2.7II2.9|) 

I • w • 02X 



r=l 



sinh(a — ~ ivr/? 
sinh(a — 



+ 



N 

+ 4 



In 2 sinh ( — a + 



+ In 2 sinh 6 + a 



-20 



. (3.18) 



Again, in the continuous limit the last term is zero in the interval < < 1/2 and we 
use our regularisation procedure. According to the latter we obtain a relation 



F-(a) = -F+ (a-m/?2) 



(3.19) 



which gives F~(a) for any complex a. From the knowledge of F^(a) it is simple to 
reconstruct A(a) = exp[F+(a;)] + exp[F^(a)] in all the a-complex plane. 

From the classical results jT2] we can argue that the coefficients of the asymptotic 
expansion of In A(a) are proportional to the local integrals of motion in involution "i^n+i 
of perturbed conformal field theories. Therefore, using the same techniques as in P], 
we expand InA(a) for Rea ±00 in the strip max| — |,— 7} < Ima < min||,7} 
(which is only a technical assumption): 



In A(a) = — mR ctg — cosh —a ± 
^ ^ 27 7 



' n=0 



T-a(2n+l' 



/+00 
'OO 



(3.20) 



Now, we need to compare this expression with the conformal analogous (5.48) of jilj and 
to consider the conformal normalisation coefficients c„ ((5.50) of P]), to write down 



In A(a) 



27 



mR ctg — cosh —a ± — ^1°'^'^^'^'^'^ Cn+ii^n+i 1 



7 



(3.21) 



n=0 



where the local integrals of motion have the exact expressions {x = -pl): 

1 TT 



■•-271+1 



c„+i 7 
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(3.22) 



Indeed, we have renormalised them with respect to the conformal counterparts hn+i 
(see (5.51) of [1 ) so that in the conformal hmit 

/ f \ 2n+l _ 

lim (— j IL"-^! = hn+i , (3.23) 

where r = mR is a dimensionless parameter. Moreover, we have been able to build up a 
lattice field theory for which conjecture (3.31) of [T3| - suggested by pure analogy with 
the conformal case [T^ - holds. Nevertheless, the naive scaling limit of our monodromy 
matrix is absolutely different from the proposed (but not actually used) prescription 
given in JH]- On the other hand, we are in agreement with the finding of [TS] derived 
using a brilliant development of the Thermodynamic Bethe Ansatz (TBA) technique 
in the case of Lee- Yang model |15j . 

For jS"^ = 1/2 the local integrals of motion (j3.22j) are explicitly given by 

/ + 00 
^^g-x2{x+i^){2n+l) In |1 + ^^[rnRsmh2{x+^r|)+2neu.] | ^ 
-oo 

(3.24) 

where < rj < 7r/2. The choice rj = tt/A gives 

iL+i = x(-l)"2^-"7r2"(n + l){2n + l)r~'"-' ReP2n+i(a, r) , (3.25) 
where a = — 2718(00 + 1/2) and we have introduced the rescaled charges 

r+co 

^2n+i(a,r) = r2"+i / dO coshe(2n + 1) In [l - e-'^"'"^^-'"] , (3.26) 
Jo 

which have a finite conformal limit (r 0). Their derivative with respect to r is 
written, after an integration by parts, in terms of modified Bessel functions. 



dP2n+l ( „ ^\ _ „2n+l „-iak 
k=l 



dr 



a, r = r 



^J2e-''''K2n{kr). (3.27) 



After integrating this formula and considering that the integration constant is fixed by 
the known conformal limit ()3.23p . we would obtain an expression of l2n+i series 
of modified Bessel functions. Nevertheless, we prefere to differentiate ()3.27|) and write 
down a relation involving Kq only : 

/ 1 J \ 2n+l +00 ,2n +°° 

-— ReP2n+i(a,r) = Y cos ak Ko{kr) = {-ly -^y cos akK^{kr) . 

\rdrJ j-^ da ^ 

^ ' fe=i fe=i 

(3.28) 

This expression allows us to simplify the form of the off-critical charges via a Schlomilch 
formula ^T] (C is the Euler-Mascheroni constant): 

00 1 r vr 

EK^ikr) cos fca = - fc + In H , ^ + 



k=\ 

00 



vr \ - 



+ 



2Vr^ 
1 1 



v/r2 + (2/7r - a)2 + {21tx + af vr/ 



(3.29) 
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Term by term derivatives can be performed using the formula 

2n 

(^2" 1 



^ b a x 



h=() 

he2N 



2 



f,(2n) 



(_l)n+t [(2n- 1)!!]' 



(3.30) 



/i! = 



We give the exphcit result for the case n = 0, 



RePi(a,r) = — [2C - 2 ln47r + Inr^ - 1] + -Vr^ + a2 - -\a\ + (3.31] 
8 2 2 



+ 



(2/7r - a)2 + v/r2 + (2/7r + a)2 - 47r/ 

27r/ 



7r25, ( — 
27r 



because in the cases uj = 0, ±1/2 it reproduces known TBA results |16j . 

So far, we have studied the monodromies of l^n+i ^ = Actually, it is also 
important to see how the charges depend on the twist u. Let us extend ()3.24|) to 
semi-integers n > 0. Introducing the index A; G N, A; > 1, we find that the following 
differential equation holds: 



V27re/ 2V27r(A; + 3) duj \2neJ2y^k{k + l)du ^+ 



0, 



(3.32) 



In the conformal limit r 0, taking into account ()3.23|) . we get the solution J^+i 
-e2-'-'/^Bk+2{uJ + l/2). 



4 Summary and outlook 

Starting from the lattice theory of two coupled (m)KdV equations, we have found in the 
whole interval < (3"^ < 1 the nonlinear integral equation which describes the vacuum 
sector and we are in the position to work out all features of excited states after j3] . The 
existence in our approach of a transfer matrix (which generates the local integrals of 
motion in involution) has allowed us to give expressions for its eigenvalues on the twisted 
vacua in terms of the solutions of the nonlinear integral equation and makes possible, in 
the next future, to find a deep link with the Thermodynamic Bethe Ansatz technique 
of jni- We have found strong evidence that local charge eigenvalues coincide with 
the primary state eigenvalues of the local integrals of motion of $1^3 perturbed minimal 
conformal models, the primary weight being linked to the twist. The expressions become 
explicit for (3'^ = 1/2. Eventually, we want to emphasise that the wide generality of 
our construction |0 |Tj makes possible its application in the domain of two-dimensional 
supersymmetric field theories. 
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